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Drops of active liquid crystal have recently shown the ability to self-propel, which was associated 
with topological defects in the orientation of active filaments [Sanchez et al , Nature 491, 431 
(2013)]. Here, we study the onset and different aspects of motility of a three-dimensional drop of 
active fluid on a planar surface. We analyse theoretically how motility is affected by orientation 
profiles with defects of various types and locations, by the shape of the drop, and by surface friction 
at the substrate. In the scope of a thin drop approximation, we derive exact expressions for the flow 
in the drop that is generated by a given orientation profile. The flow has a natural decomposition 
into terms that depend entirely on the geometrical properties of the orientation profile, i.e. its bend 
and splay, and a term coupling the orientation to the shape of the drop. We find that asymmetric 
splay or bend generates a directed bulk flow and enables the drop to move, with maximal speeds 
achieved when the splay or bend is induced by a topological defect in the interior of the drop. In 
motile drops the direction and speed of self-propulsion is controlled by friction at the substrate. 

PACS numbers: 87.16.-b., 82.70.-y., 68.15.+e 


I. INTRODUCTION 

Active fluids are non-equilibrium materials composed 
of elongated units that constantly consume energy to 
move themselves forward or to stir their surrounding 
medium m- Living systems that can be described as 
active fluids are often found in spatial confinement, like 
suspensions of swimming bacteria in porous media [4] or 
the cytoskeleton [5 of a cell enclosed by its membrane. 
Experiments with active fluids confined to droplets re¬ 
vealed a myriad of novel features, such as spontaneous 
symmetry breaking |6], self-organised defect structures 
mm®, coherent large-scale flows M2\ and even self- 
sustained motility [T . A full theoretical understanding 
of the interplay of activity with total geometrical con¬ 
finement is still lacking, yet it could not only help to un¬ 
veil the basic fluid dynamical processes underlying cell 
motility and cell division, but also allow to tune the 
macroscopic behaviour of biomimetic droplets [lj 8], for 
instance by controlling the orientation of active filaments. 

Previous theoretical work on active droplets has 
focused on two-dimensional geometries or spherical 
droplets immersed in an other fluid. The active fluids 
that these models typically refer to are aqueous suspen¬ 
sions of cytoskeletal filaments together with their asso¬ 
ciated motors, in particular actin filaments with myosin 
motors as a material having polar symmetry or micro¬ 
tubules with kinesin motor clusters [T] as one having 
nematic symmetry, respectively. In a circular domain 
an active fluid can self-organise into a stable circulat¬ 
ing state at high enough activity m, which is related 
to cytoplasmic streaming in plant cells. A similar rotat¬ 
ing spiral, albeit with a counterrotating boundary layer 
and varying cell orientation, was observed in a quasi- 
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two-dimensional circular suspension of swimming bacte¬ 
ria [9]. Hydrodynamic interactions among the bacteria 
and with the confining wall were identified to be crucial 
to drive this self-organised pattern [10]. Nematic two- 
dimensional droplets showed elongation, spontaneous di¬ 
vision and motility as a result of the interplay of defects, 
geometry and activity El- Finally, in a numerical study 
of polar active droplets of initially circular (2D) or spher¬ 
ical (3D) shape, surrounded by a passive fluid, increasing 
activity led to spontaneous symmetry breaking accompa¬ 
nied by self-propulsion and deformation p~5l [16] . These 
findings were underpinned theoretically by the analysis 
of a droplet of fixed shape and with an imposed splayed 
orientation field m A related class of models is con¬ 
cerned with the cellular actin cortex as an active gel and 
provided insight into activity-driven deformation and mi¬ 
gration of cells HBHH]. 

The behaviour of a three-dimensional active drop 
placed on a fiat surface l2D423j has been hitherto scarcely 
investigated. However, it is important to understand 
this setup in order to clarify the role of active flows in 
motile cells on substrates [24] and, in a broader fluid dy¬ 
namical perpective, to characterise novel phenomena not 
found in the well-studied passive fluid droplets on sur¬ 
faces [25]. Topological defects in the orientational order, 
which are points of undefined orientation [26] . emerge 
spontaneously in active systems [lj. Due to strong gradi¬ 
ents of active stresses in their vicinity, defects can act as 
sources for large-scale flows which, in combination with 
the confined shape, can drive macroscopic motion of a 
drop of active fluid along a surface. To analyse which 
types of orientation profiles with defects lead to three- 
dimensional flows that drive a drop from within, and how 
surface friction at the substrate and the shape of the drop 
affect these flows, is the purpose of this work. 

A first theoretical analysis of active drops on surfaces 
focused on spreading laws and stationary shapes for sev- 
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eral polarisation fields with high symmetry m , account¬ 
ing for a no-slip boundary only. In a recent numerical 
study of a drop of polar active fluid various biologically 
relevant shapes and motile steady states were obtained 
[22] . The computation included both actin treadmilling 
and a variable effective friction representing the focal ad¬ 
hesions as mechanisms that are crucial for cell motility 
on a substrate. In another computational study realistic 
three-dimensional shapes of motile cells were generated 
using a hydrodynamic model [23] . although without in¬ 
cluding the polarisation field of active filaments. 

In this paper we consider generic flow profiles that 
arise in three-dimensional droplets on planar substrates 
with a variety of polarisation fields associated to differ¬ 
ent types of defects. We derive analytical expressions for 
the flow field resulting from a specified polarisation and 
varying friction at the substrate. This leads to a geo¬ 
metric description of the flow, with a natural separation 
of splay and bend into horizontal and vertical compo¬ 
nents, of which only the latter depends on the shape of 
the drop. We find that self-propulsion of the drop is en¬ 
abled by asymmetric orientation fields and is fastest when 
the splay or bend is induced by a defect approximately 
mid-way between the centre and boundary of the drop, 
for instance an aster or a vortex. When the substrate 
friction is negligible the interior flow in the drop is rota¬ 
tional and there is no propulsion. Maximal propulsion is 
achieved when the friction is greatest. The shape of the 
drop is free to deform as a result of the active flows. We 
qualitatively describe the deviation from an initial spher¬ 
ical cap shape for two examples. For the case of an aster 
defect in the centre of an axisymmetric drop we derive 
the stationary shape of an active drop. 

The paper is structured as follows. In Sec. [!!]the model 
for a thin drop of acti ve flu id is derived, including the 
equations for the flow (II A), the most general for m of a 
polarisation field subject to tangentia l anc horing (IIB), 
and the appro ximat ion of equations (IIC) and bound¬ 
ary conditions (IID) under the assumption of separable 
length scales. The exact general solution for the flow 
field is presented at the start of Sec. m In Sec. m A] the 
structure of the horizontal flow is explained and illus¬ 
trated by examples of various orientation fields. In Sec. 
iBBl the implications of the surface friction on motility 


are discussed and Sec. Ill C| is concerned with shape de¬ 
formations and stationary shapes. Finally, a discussion of 
the results and their relation to experiments is provided 
in Section HVl 


II. MODEL 

We consider a three-dimensional drop of active fluid, 
illustrated in Fig. [l] with mixed boundary conditions: a 
flat, rigid substrate providing frictional dissipation un¬ 
derneath and air or a fluid of much lower viscosity with 
zero tangential stress above the drop, allowing the shape 
to deform. The active fluid in our model is either an ac¬ 



FIG. 1. (a): Schematic view of a droplet of active fluid with 

filaments oriented according to the field S', which has the 
two-dimensional projection P. The filaments are anchored 
parallel to the bounding surface, which is given by its height 
function h. (b): The active stresses induced by the filaments 
generate a flow u. 


tive liquid crystal [27] or an active polar gel [28]. repre¬ 
senting a suspension of active filaments in nematic or po¬ 
lar long-range orientational order, respectively. We study 
how the onset and different features of motility of such 
drops are controlled by spatial variations and defects in 
the orientation field, the shape of the drop, and the sur¬ 
face friction at the subtrate. To address these questions 
we derive a simple hydrodynamic model, which is built 
on two basic assumptions. 

Active fluids are described theoretically in terms of 
generalised hydrodynamic equations that couple the dy¬ 
namics of the flow and the orientation field and can be 
derived by coarse-graining a suspension of active parti¬ 
cles mmm or by modelling the cytoskeleton in terms 
of a viscoelastic gel driven out of equilibrium [281 EH [32] . 
The equations are a coupled non-linear system for both 
the orientation field of the active filaments and the ve¬ 
locity field of the active fluid. A simplification can be 
obtained for dynamic steady states by assuming that the 
alignment is static and close to the energy-minimising 
configurations of passive liquid crystals. That given, one 
need only solve for the active fluid flow. We impose tan¬ 
gential anchoring at both the base of the drop and the free 
surface, while at the contact line the orientation is free 
to point in any direction. This choice of anchoring is mo¬ 
tivated among others by experimental setups where the 
filaments adhere to the boundary layer in large enough 
droplets BE]. Also, in the lamellipodia of crawling cells 


















3 


actin filaments lie tangential to the substrate and push 
perpendicularly against the leading edge [33] . 

The second assumption is that the average height of 
the drop is much smaller than its width and we can sepa¬ 
rate the two length scales to simplify the equations for the 
flow field, as is typically done for thin fluid films mm- 
For this purpose we choose the drop to be initially in the 
shape of a flat spherical cap with radius and height 
function given by z = h(x,y) = y/l — x 2 — y 2 — 0.5 (see 
Fig. [lja)). However, the shape is free to deform due to 
the active flows generated in the drop and these defor¬ 
mations will also change the orientation field indirectly, 
through its anchoring to the surface. 

These assumptions allow us to find an exact solution 
for the flow field in a drop of a given shape and with 
a prescribed orientation profile. Its structure turns out 
to be guided by both the familiar instabilities towards 
bend and splay in active systems m and the coupling of 
the orientation to gradients in the height profile, which 
has previously been studied in the context of spreading 
of drops [21]. The exact formula for the flow allows us 
to address a variety of aspects of drop motility in an 
analytical way. Our model reveals that asymmetric spa¬ 
tial variations in the orientation are key to enable drop 
propulsion and explains why friction at the substrate is 
essential. We also find that topological defects can both 
drive drop motility and also induce singular vertical flows 
that could lead to the growth of protrusions at the site 
of the defect. 


A. Governing equations 


Here, we only consider the dynamics of the flow field 
u{x , t ) = ( u , v, w) and the pressure p(x, t ) of the suspen¬ 
sion, which comprises both the active particles and the 
solvent fluid as a whole, and the resulting variations in 
the height profile of the drop, given by z = h(x,y,t). 
The hydrodynamic effect of the active filaments on the 
suspension is an active stress, whose gradient drives the 
flow. To first order in a gradient expansion the active 
stress tensor for nematic and polar suspensions is indis¬ 
tinguishable 21. EZ1, which is applicable here since we 
are interested in flows that are of large scale compared 
to the size of a filament. As the systems described here 
are in the low Reynolds number regime, the dynamics 
of the suspension is given by the generalised Stokes and 
continuity equations [21] 

- Vp + /jlA u + V • (cr a + cr e ) = 0, (1) 

V-u = 0, (2) 


where p is the density, p the viscosity, cr e the Ericksen 
stress of passive liquid crystals, and cr a the active stress. 
The boundary conditions for equations 0 and 


are 


specified below, in Sec. IIP 


B. Polarisation field and active stress 

The polarisation is taken to be tangential at both 
bounding surfaces, the flat base in contact with the sub¬ 
strate and the free upper surface of the droplet, z = 
h(x,y,t). If the polarisation on the base (z = 0) is 
(P x ,P y ,0), with \P\ = 1, then this can be matched onto 
a tangential polarisation on the upper surface by setting 

Pz = (PAh + Pydyh)/^ 1 + {d x Kf + (d y /l) 2 . (3) 

The resulting vector field lies tangential to the 
surface with the unit outward no rmal rih = 
(■ -d x h , -d y h, l) T /yi + (d^/i) 2 + (dyhjz. Finally, the 
three-dimensional polarisation field S(x) is constructed 
by smoothly interpolating between the two surfaces using 
the factor z/h in the ^-component, 

i ( hP x \ 

s{x> - sw et j • <4> 

and normalising with n(x) = %/h 2 + z 2 P' 2 . This is the 
most general form of a polarisation field that satisfies the 
tangential anchoring condition and is consistent with the 
assumption of a thin drop, with small variations of the 
polarisation in the ^-direction. 

The active stress tensor is defined as m 

°ij = -00 (.SVS'j - y) » (5) 

where the constant <ro sets the strength of the activity 
and whether the stress is contractile (<to < 0) or extensile 
(ao > 0), that is whether the active particles drive fluid 
in or out along their axis of orientation, respectively. 


C. Separation of length scales 

We can exploit the geometry of the drop to simplify 
equation 0- We assume that its horizontal extension, 
characterised by the length scale L, is much larger than 
its average height ho. Therefore, as we are interested in 
long-wavelength phenomena, the variations of the flow 
in x and y directions are much more gradual than in z 
direction. We define the small parameter s = ^ 1 and 

perform an expansion of the governing equations, similar 
to the approach for thin films [2H [25U34]. 

We make the coordinates dimensionless by scaling ac¬ 
cording to x —y 7/ —z and likewise h -A . 

Compared to the°height tiie extension of the drop is sim¬ 
ilar in both horizontal directions, therefore we treat the 
x and the y directions analogously. With a characteristic 
velocity Uo we scale u ^ and v jj~. For the ver¬ 
tical flow component the incompressibility condition ([2| 
requires w —>> Finally, the time scale of the flow is 

set by L/Uq, and therefore t —)> 
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When we expand equation 0 in 5 then, in the x com¬ 
ponent for instance, the leading order term is ^ d\u , 
as expected for a thin film [25] . while the contributions 
due to pressure and active stress are of sub-leading order. 
Therefore, to retain the effect of activity these terms are 
scaled as 


sh 0 eh 0 

V -> and cr 0 ->► -^ctq. 

yU 0 vUq 


(6) 


The field P is normalised, so its components are P Xl P y ^ 
0(1). It follows that the ^-component 0 should be 
rescaled as P z —>■ Pz/e and the normalisation of S be¬ 
comes n{x) ~ hoh(x,y). With these scalings the dom¬ 
inant terms of the gradient of the active stress can be 
written as (V • cr a )_\_ = ^in terms of an effective 
active force 


fl = -(70 (p( Vj_- P + • Vxft) + (P-Vl)P ) , (7) 


which only depends on the horizontal position x± = 
(x, y). We find that has a natural decomposition into 
terms associated with the geometric properties of P, that 
is its splay and bend, and a term that couples P to gra¬ 
dients in the height profile of the drop. This structure is 
inherited by the resulting flow field, which is driven by 


and will be discussed in detail in Sec. Ill A 
The tensor cr e in 0 accounts for standard nematic 
elasticity [26] . Assuming one elastic constant IT, the 
leading order term in the elastic distortion free energy 
density is the one associated with splay, y (V • S) 2 ~ 
0(e 2 ). Thus, in relation to the active stresses we can 
omit cr e . In addition, the orientation profiles we consider 
differ from the ones that minimise elastic energy for pas¬ 
sive liquid crystals only by a small vertical tilt, suggesting 
that flows due to elasticity should be small compared to 
active flows. 

To leading order in 5 equation 0 becomes 


d 2 z u i — V . p 
0 = d z p , 


(8) 

(9) 


where u± = (it, v) and Vx = {d x ,d y ), while the continu¬ 
ity equation remains unchanged. 


D. Boundary conditions 


In this section we first state the boundary conditions 
in terms of unsealed variables, before transforming to di¬ 
mensionless variables and expanding in the small param¬ 
eter 5. At the base we allow for partial slip with linear 
friction, such that the shear stress acting on the fluid is 
proportional to the fluid velocity at the boundary, 


rig Te* = 


— -iLi for i = x,y 


( 10 ) 


Here £ is the inverse friction coefficient, so that £ = 0 cor¬ 
responds to the no-slip condition. The full stress tensor 
of the fluid is given by 


-pSij + y (djUi + diUj) + afj + < 


2J5 


( 11 ) 


and no = — e z is the outward normal to the lower bound¬ 
ing surface. After making the inverse friction coefficient 
dimensionless by scaling it as £ for small e condi¬ 

tion (10) becomes 


u±(z = 0 ) 


-£<9 z u± 


z =0 


( 12 ) 


now written in dimensionless quantities. The substrate 
is non-permeable and therefore the ^-component of the 
flow must vanish there, w \ z= o = 0 . 

At the free boundary the stress component tangent to 
the surface should vanish, 


njT - (nlTn h )n h = 0, (13) 

which for small 6 and in dimensionless variables simplifies 
to 


d z u ± \ z=h = 0. (14) 

The normal stress at the upper boundary is given by 
the mean curvature of the surface V • n^, the surface 
tension 7 and the pressure po in the surrounding medium, 

P = P0 + Y~ £3 ~^ 0 ( d x h + d y h ) + O ( e ‘ 2 )^ (I 5 ) 


which in the absense of activity reduces to the Young- 
Laplace equation. In (15) only 7 still carries dimensions. 
The surface tension term is three orders of magnitude 
smaller than the leading term. If we wanted to retain 
surface tension effects, the coefficient had to scale like 
7 ~ £ -3 . However, earlier we identified that the active 
stress coefficient scales as ao ~ £ _1 (cf. eq. §)• Because 
here we are interested in active effects primarily, we will 
neglect the surface tension term. Then the boundary 
condition (15) becomes p\ z =h = Po + yp. 

The upper bounding surface is free to deform and its 
time evolution is given by the flow component normal to 
the surface 


d t h = u • rih = w — ud x h — vd y h . (16) 

This kinematic boundary condition can be expressed in 
terms of u and v only, if we integrate the continuity equa¬ 
tion (with re(0) = 0) and use Reynolds’ theorem. The 
height profile then evolves according to 
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FIG. 2. (a)-(d): Top view of drops with a splayed polarisation field (a = 0, b = 1) due to a defect at ro = —2r^, —r^, — |r< 2 , 0. 

(e)-(h): Plot of the resulting flow field (18) at z = 0.01, for a contractile drop with no slip (cro = —1, £ =_0). The flow is 
aligned with the polarisation 


but changes sign along a line where the splay of P balances the coupling term 
the following plots red lines represent the polarisation field, and the flow field is decomposed into direction (white arrows 
magnitude (colour coded, scaled). At a defect the flow typically diverges and is cut off in the display (black region) 
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FIG. 3. (a)-(d): Top view of drops with a bent polarisation field (a = = 1) due to a defect at ro = —2r^, — r*d, — |rd, 0. 

(e)-(h): Plot of the resulting flow field (18) at z = 0.01, for an extensile drop with no slip (cro = 1, £ = 0). The flow is 
perpedicular to the polarisation in the bulk, but in (e)-(g) it alignes with the polarisation close to the boundary due to the 
coupling term (21). In the axisymmetric case, (d) and (h), the coupling term vanishes and the flow points radially outwards. 
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III. RESULTS 


We solve the equations <© and <© for the instanta¬ 
neous flow field u resulting from a particular orientation 
field P and shape h. From © we find, that the pres¬ 
sure is constant to first approximation, p = p 0 + 
whereby the surface tension was neglected. Then, in¬ 
tegrating equation (|8| twice gives an expression for the 
horizontal flow components 

(z 2 

u± = a ° v T + h ^ ^ 

x(p(Vx- P + tp • Vx/fj + (P-Vx)P^ (18) 


The ^-dependence of u± is primarily determined by the 
effective active force /J ( see eq. ([7])), and can be decom¬ 
posed into three parts, each with a clear interpr etatio n 
in terms of spatial variations in P and h (see Sec. Ill A). 
The vertical flow component follows from the continuity 
equation ©> 


W = -gVx-/i- 


-6* Vx-(/i/l) 


(19) 


from which it is obvious that the three-dimensional na¬ 
ture of the flow will be most prominent in regions of 
strongly varying for instance at topological defects 
in P. Fig. [ljb) provides an insight into one example of 
the full three-dimensional flow field (u±,w), however in 
the following we restrict the visual presentation to two- 
dimensional cross-sections. 

We illustrate the properties of the flow in the drop 
on examples of polarisation fields that are generated by 
defects of different topological strength and are given in 
the form 

P = (cos (a + bO(r 0 )), sin (a + bO(r 0 ))) , (20) 

where b is the strength of the defect and a controls the 
shape of the field around it, for example whether it is 
an aster or a vortex for strength +1 (see Fig. [2|d) and 

Fig. |s|d)). The angle 0(ro) = arctan (^ x l rQ ^ depends on 

the position ro of the defect along the x-axis, such that 
the defect can be located inside or outside of the drop 
(see Fig. [2|a) and (c)). We remind the reader that all 
results are presented in terms of dimensionless quantities 
introduced in Sec. 


A. Splay and bend drive directed large-scale flows 
in the drop 

Writing the effective active force, and thus the flow 
field, in a coordinate free form (see eq. @) allows us to 


interpret the contributions to the flow in terms of the ge¬ 
ometrical properties of P and h. The first P-dependent 
term in ( p~8| ) produces flow in direction of the polarisation 
field and is proportional to the sum of its splay, Vj_ • P, 
and the term 


U h = —P • \7±h 
h 


( 21 ) 


which couples the polarisation of the filaments to the 
shape of the drop and accounts for the splay into the 
third dimension. Since term (21) is large in the vicinity 
of the contact line, where h 0, it dominates the flow 
at the boundary and aligns it there with the polarisation 
field. The second contribution, (P • V_l)P, creates flow 
that is perpendicular to the polarisation field and equal 
to its bend. This is a manifestation of the instabilities 
towards splay or bend in contractile or extensile active 
fluids, respectively [3 . For a drop this mechanism means 
that both pure splay or pure bend, in an asymmetric 
configuration, can generate a directed flow in the bulk 
of the drop and thus enable it to propel itself along the 
substrate. The direction of propulsion depends on the 
sign of the activity ao and on the friction parameter £ (see 
Sec. Ill B| ) . Figures [2] and [3] show examples for orientation 
fields with pure splay or pure bend and the resulting flow 
fields. 

A polarisation field with pure splay can be produced 
by varying the position of an aster defect, as illustrated 
in Fig. [2j The resulting flow will be aligned with the 
polarisation and have a large component in direction of 
the splay. However, the flow changes its direction along 
a line where the splay of P is balanced by the coupling 
term (21), Uh + V_l • P — 0. At the right boundary the 


flow is driven by the vertical splay due to the tangential 
anchoring to the bounding surface and directed inwards 
in this example. Note, that if P is uniform then the 
coupling to the surface, that is term (21), is solely re¬ 


sponsible for generating a flow in the drop, which then 
leads to symmetric spreading studied extensively in m- 
Here however, the flow at the boundary counteracts the 
bulk flow that is driven by horizontal variations in P. 

A polarisation field with pure bend, as in Fig. [3J 
can be generated by a variably positioned vortex defect. 
The flow is perpedicular to the polarisation in the bulk, 
having a large component in direction perpendicular to 
the bend, but alignes with the polarisation close to the 
boundary due to the coupling term (21). However, in the 
axisymmetric case in Fig. j3^h) Uh vanishes, because h is 
now constant along closed circles in P, leading to purely 
radial flow. In droplets of actin and myosin solution a ra¬ 
dially inward flow, that was associated with an emergent 
ring-like structure of actin, was observed experimentally 
nu- We can recover this result by choosing contractile 
activity ( ctq = — 1) in Fig. |3jh). 

For more complex polarisation fields all contributions 
to fll8| ) are present, as shown in Fig. [ 4 ] for three different 
types of defects placed in the centre of the drop. The 
spiral defect (Fig. Qa)) generates rotational flow around 
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FIG. 4. Left: Top view of polarisation fields with (a) spiral 
(a = 71 / 4 , b = 1 ), (c) + 1/2 defect (a = 0 , b — |), and (e) — 1/2 
defect (a = 0,6 = — |) in the centre of th e dr op (ro = 0 ). 
Right: Plot of the corresponding flow fields |l 8 | for z — 0.01. 
Parameters are ao = 1 (extensile) and £ = 0 (no slip). 


the centre. However, close to the boundary the direction 
of rotation reverses, which is here again a result of the 
coupling to Vift. A similar flow pattern, including the 
change of rotational direction, was observed in thin drops 
of bacteria suspensions where the swimming bacteria self- 
organised into a spiral vortex [9|. 

In addition to profiles with integer strength defects, 
which can be found in systems with polar or nematic or¬ 
der, we also consider half-strength defects that only exist 
in nematic systems. Such defects emerge spontaneously 
and shape the dynamic flows in many active systems with 
nematic symmetry mm, particularly the self-propelling 
+1/2 defect. In our model, a +1/2 defect in the cen¬ 
tre of the drop (Fig. (4^c)) creates unidirectional flow in 
the bulk, unlike +1 defects in the same position, because 
it has an advantageous combination of a splayed and a 
bent region that drive the flow in the same direction. A 


— 1/2 defect (Fig. [4je)) on the other hand leads to no 
net direction of the flow in the bulk due to the three-fold 
symmetry of the corresponding polarisation field. These 
flows are in agreement with the two-dimensional results 
for half defects, that explain the self-propulsion of posi¬ 
tive and the stagnation of negative half defects in active 
nematics [35, 36]. 

A point defect in the field P, when positioned in the 
interior of the drop, corresponds to a line defect in the 
field S. At this line defect the magnitude of the horizon¬ 
tal flow typically diverges, for instance as ^ ^ at the two 
+1 strength defects in Fig. [2|d) and[3]]d). This could re¬ 
sult from the fact that in the very vicinity of a defect the 
approximation of separable length scales breaks down, 
because the flow varies strongly on a small length scale 
in the horizontal direction, as well as in the vertical, and 
all spatial derivatives should remain in eq. 0 . However, 
our model offers a realistic prediction for the direction 
of the three-dimensional flow at defect lines. In particu¬ 
lar, the radial inflow at +1 defects is converted into the 
third dimension by continuity and could lead to a thin, 
protrusion-like deformation of the drop at the location 
of the defect. Such protrusions were recently observed in 
experiments of thin shells of active liquid crystal [ 8 ], un¬ 
der the condition of lowered surface tension with respect 
to the filament elasticity, and could be the experimental 
equivalent of such singularities in flow that are induced 
by a defect. 

So far we have excluded the effect of surface tension 
on the flow in the drop, because the associated term was 
of lower order in the approximation. If one chooses to 


include this effect and scale 7 ~ £ 3 in eq. (15), it will 
contribute 


U surf = - ( 7 + h g - z) ) 7 V±Vi h ( 22 ) 


to the horizontal component of the flow. For a no-slip 
boundary, this would yield an additional radially inward 
flow and thus either reduce or enhance the flow at the 
boundary due to ©, depending on the sign of ctq. Thus, 
the bulk flow would still be primarily determined by vari¬ 
ations in the horizontal orientation profile. 


B. Surface friction controls speed and direction of 
motile drops 

To determine whether a drop, which produces a di¬ 
rected horizontal flow in the bulk (like in Fig. [2jf ) or 
|3]]f)), will move its centre of mass it is necessary to in¬ 
vestigate the vertical flow component. The amount of 
slip at the rigid surface, represented here by the effec¬ 
tive friction parameter £, controls whether the flow in 
the bulk of the drop is laminar or rotational, and in the 
former case determines the direction of the flow. 

Thus, surface friction has implications for the direction 
and speed of the centre-of-mass movement of the drop, 
which is summarised in Fig. [5] for the example of a splayed 

























FIG. 5. Centre-of-mass velocity for different values of fric¬ 
tion and different positions of an aster defect in a contractile 
drop (or of a vortex defect in an extensile drop, respectively), 
with marked maxima (red dots) and minima (blue squares). 
A no-slip boundary yields the highest maximum velocity in 
the direction away from the defect and for large enough slip, 
£ > 0.2, the direction of propulsion is reversed. Maximal 
propulsion speeds are achieved with the defect being placed 
asymmetrically in the interior of the drop, — 1 < ro < 0. 


orientation field in a contractile drop. As a measure for 
the self-propulsion of the drop along the substrate we 
numerically calculate the centre of mass velocity v cm as 
the integral of the flow component in x-direction, v cm = 
Fo idrop^F. This velocity depends on the friction and 
the defect position, and reaches a maximum when the 
defect is located inside the drop. Exactly the same plot 
holds for a bent orientation field in an extensile drop. 
Surface friction is essential for the self-propulsion, since 
a no-slip boundary yields the highest possible maximum 
velocity, while an increasing slip brings the drop to a halt 
and, for £ > 0.2, reverses the direction of propulsion. 

We illustrate the effect of friction on the flow by the 
example of the splayed orientation field from Fig. [2jb) , 
where an aster defect is located on the contact line. For 
a no-slip boundary (£ = 0) the flow vanishes at the base 
and, at finite height, is directed away from the defect 
in a mostly laminar way in the bulk of the drop (not 
shown). For small slip (Fig. [6ja)) a thin treadmilling 
layer emerges close to the base, on top of which the flow is 
still laminar. Friction with the substrate induces a shear 
flow that opposes the flow generated by active stresses. 
For finite slip this creates a vortex in the fluid that moves 
upwards and spans a larger region of the drop with in¬ 
creasing slip (Fig. [6^b)). If the friction is small enough, 
the flow becomes laminar again but with reversed direc¬ 
tion compared to the no-slip case (Fig. §<=))• 

This behaviour is apparent from the role of the friction 
parameter in the ^-dependent prefactor in ( p~8] ) , which is 
bounded from below and above by 


h U - < y + h (£ - z ) < • (23) 


While the upper bound is h £ > 0, the lower bound be¬ 



z 0.2 


FIG. 6. Plot of the flow components (u, w) for y — 0, given by 
(18) and (19), for the splayed polarisation field in Fig. [2^b) in 
a contractile drop (cro = — 1). In this side view the defect is 
located in the left corner. The slip increases from (a) to (c), 
£ = 0.07,0.15,1.5. The color-coded magnitude of the flow 
corresponds to (c) and scales in the same way as £. 


comes negative in those regions of the drop where 

t h 
2 • 


(24) 


Both horizontal flow components u and v then change 
sign at a height zo(x,y) = h — y/h(h — 2£). For high 
enough friction condition ( [24] ) is satisfied in the bulk of 
the drop and the flow is rotational. In the low friction 
regime the transition from rolling to spreading upon in¬ 
creasing slip is analogous to what is found for passive 
fluid drops on an inclined surface m, where a longer 
slip length corresponds to more sliding and less rolling of 
the drop. 

Fig. § a) represents a realistic scenario for the bulk 
of a polarised and motile cell extract of actomyosin en¬ 
closed by a membrane [38] , where the effective friction is 
mediated by focal adhesions that attach the actomyosin 
gel and membrane to the substrate. Here the bulk of 
active fluid flows in direction of the splay, away from 
the defect, thus such a cell extract would move to the 
right. The backslip of membrane at the substrate could 
be accounted for by a reduced number of focal adhesions. 
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FIG. 7. Stationary shape profile h(r) for an axisymmetric 
drop of unit volume with an aster defect in the centre (see 
Fig. |2] (d)) and three different final heights ho < 3£, £ = 1. 


This flow field is consistent with the forward flow of cy¬ 
tosol observed experimentally in a moving cell viewed 
from above [24 and with with numerical results for the 
three-dimensional flow in a crawling cell [ 22 ] . 

The backward flow at the leading edge in Fig. | 6 ][a) 
(see Fig. If) for top view) is a result of the strong tan¬ 
gential anchoring of filaments to both bounding surfaces, 
which splays the filaments vertically. This flow enhances 
spreading of fore-aft symmetric active droplets [ 21 ], but 
in the case of a strong directed bulk flow driven by hor¬ 
izontal variations in the orientation, as in Fig. [ 6 ] it can 
turn into a backflow on one side and counteract drop 
propulsion. There are two ways in which our model 
could be modified to eliminate this backflow at the lead¬ 
ing edge. Firstly, we could locally remove the tangen¬ 
tial anchoring condition in a region opposite the defect 
and with it the source of the backflow. Secondly, we ex¬ 
pect that including self-propulsion of the filaments along 
their direction of orientation, which is the simplest way 
to model actin treadmilling Hang, would enhance the 
bulk flow and compensate the inward flow at the frontal 
boundary. 


C. Shape deformations 


The condition of free tangential stress at the upper 
bounding surface means that this boundary is free to de¬ 
form in response to the active flows pushing against it. 
The changes in shape can be expressed in terms of the 
horizontal flow components (see eq. ([T7|)) , which them¬ 
selves depend on the shape. The resulting self-consistent 
equation for the height function allows to not only find 
the immediate shape deformations due to a particular 
orientation field P, but also to look for steady state so¬ 
lutions where the shape is adjusted such that the given 
P does not induce any deformations of the drop. 

In the axisymmetric case of an aster defect line in the 
centre of the drop (Fig. [2|d)) the time evolution of the 



FIG. 8. Shape deformation away from the spherical cap in 
the units of a small time step At for a contractile drop with 
(a) splayed and (b) bent orientation field (see Fig. [2ja) and 
Fig. [3][a)) and a no-slip boundary. 


shape is given by 


dth = cro - d r 
r 


£ - 3 ) hd r ( rh ) 


(25) 


For finite slip, £ > 0, this can be solved to give an ex¬ 
act, implicit solution for the stationary shape of a drop 
of radius and height ho at the centre. The constraint 
of finite volume is approximated here by the condition 
hor d = 1. The radial dependence of the axisymmetric 
steady state profile is shown in Fig. [7] for small heights 
ft-o < 3£. In particular for small values of ho our solution 
is very similar to the “flat pancake” shape that was ob¬ 
tained numerically in m for small pressure. The contact 
angle 0 for drops with ho < 3£ is given by 


tan 0 = 


_ 2/tg(3g - hp) _ 

h 0 + 3£ + y3(3 £ 2 + MX - ho)) ' 


(26) 


It takes the maximum value of # max = 7 t/4 for ho = 
2 £, which also marks the cross-over from rotational to 
laminar flow in the large slip regime (cf. eq. (|24|)). For 
ft-o > 2 £ the shape function h(r) develops a singularity 
at the origin and for ho > 3£ the solution we get is not 
consistent with a finite drop, indicating a lower limit for 
the amount of friction required to obtain a stationary 
shape profile. 

We note that the steady state shape for a vortex defect 
line (Fig. |3jd)) in the centre of the drop is given in | 2 T] as 
a non-mo notonic function with a minimum at the origin. 
We obtain a flow field that is consistent with this shape 
for an extensile suspension with a no-slip boundary, if we 
include the surface tension (see eq.(|22|)). In this case the 
defect generates radially outward flow, which is opposed 
by an inward flow at the boundary due to surface tension, 
suggesting that the fluid will accumulate along a ring. 

For a general polarisation field P the drop shape 
evolves according to 


d t h = crodi 



3 ) hd 3 (hPiPj) 


(27) 
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For the asymmetric polarisation fields considered in III A 


it is not possible to obtain an analytical result for the 
steady state. We can however qualitatively describe 
the deviation Ah from the initial spherical cap shape 
at time to after a small time interval At if we plot 
Ah « dth(to)At for a particular orientation field (Fig. 
[8|. We consider a contractile drop a with no-slip bound¬ 
ary for this purpose. For a splayed orientation field as 
in Fig. |2ja) the instantaneous flow results in an increase 
in height along a crescent-like area and a decrease to the 
front and back of it (Fig. [8^a)). For a bent orientation 
field as in Fig. [3ja) the height grows in the rear of the 
drop and descreases along the sides, so we can expect the 
drop to become more slender. However, the changes in 
shape will alter the orientation field and thus the flow. 
To investigate the dynamics of the shape over longer time 
periods it will be necessary to perform numerical simu¬ 
lations, which is not in the scope of this paper. 


IV. DISCUSSION 


We have given an analytic description of dynamic 
steady state motion of active drops. In the scope of a 
thin drop approximation, we derived exact expressions 
for the flow in the drop that is generated by a given ori¬ 
entation profile. This leads to the identification of two 
key requirements for the self-propulsion of thin droplets 
of active fluid on planar substrates. These are an asym¬ 
metrically splayed or bent orientation field, for instance 
induced by a topological defect in the interior of the drop, 
and sufficient surface friction provided by the substrate 
for the drop to push itself forward. The active flows in 


the drop are primarily driven by horizontal variations 
in the orientation and the coupling to the shape con¬ 
tributes mainly at the contact line, where it can generate 
flows that work against drop propulsion. A motile drop 
will become stagnant, treadmilling on the spot, when the 
surface friction drops below a certain value. A substrate 
with a large slip length could produce drops that move 
in a direction opposite to the expected, which is in the 
direction of splay for contractile cell fragments [38] or in 
direction of the comet-head for drops of extensile active 
nematics with a positive half-defect [T]. 

Our model suggests that the flow of cytosole in a motile 
cell fragment can result from active stresses in the po¬ 
larised and splayed actin cytoskeleton. The result is 
in qualitative agreement with the forward directed flow 
measured experimentally in rapidly moving keratocytes 
[24] . Although the role of fluid flow in the crawling ma¬ 
chinery is still unclear [39] . it was suggested to be re¬ 
sponsible for active transport of actin monomers to the 
leading edge uni. 

Shape deformations of active droplets can be studied 
experimentally on vesicles filled with an active suspen¬ 
sion, where the surface tension can be reduced by apply¬ 
ing hypertonic stress [8]. To impose a particular, static 
orientation of active filaments the drop could be placed 
on a micropatterned substrate m, specifically to induce 
axisymmetric defect configurations which would allow to 
observe the corresponding stationary shapes of the drop. 
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